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Abstract 

We propose a new approach to the Rayleigh-Schrodinger perturbation expansions of 
bound states in quantum mechanics. We are inspired by the enormous flexibility of 
solvable interactions with several (N) discontinuities. Their standard matching solu- 
tion is modified and transferred in perturbation regime. We employ the global renor- 
malization freedom of the local wave functions and derive a compact iV— dimensional 
matrix formula for corrections. In applications, our recipe is shown non-numerical 
for all polynomial perturbations of any piece-wise constant zero order potential. 
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1 Introduction 



Rayleigh-Schrodinger perturbation theory |1J leads to several popular and efficient 
numerical approximation methods [0 . Its construction of observables may also prove 
inspiring in the more abstract analysis of their coupling dependence The latter 
role of the perturbative power series ansatzs has already been emphasized in the 
classical monograph by Morse and Feshbach [Q . Their presentation of the Rayleigh- 
Schrodinger formalism contemplates any potential V(r) = V^°\r) + XV^'(r) as 
defined on a trivial square- well background V^°\r). 

In the late sixties the "mainstream" attention has been shifted to analytic V^(r). 
People have noticed that a suitable normalization leads to five-term recurrences and 
to an enormous simplification of the construction of the anharmonic oscillators with 
(r) ~ r 4 || etc. The history has been reviewed, e.g., by Simon |J. 

An unpleasant obstacle to a broader applicability of perturbative solutions formed 
by the power series in A is definitely the narrow variability of the available analytic 
zero order approximations. In three dimensions Newton |7j only lists square well 
y(°)(r) ~ (r/L) p , p -> oo, harmonic well 7 (0) (r) ~ (r/L) 2 , Coulomb field V (0) (r) ~ 
(r/L) _1 and a rather exotic strongly sing ular V^(r) ~ (r/L)" 4 §. 

The non-analytic square well seems tedious in comparison but in the present paper 
we still return to discontinuous V^°'(r). A new version of the Rayleigh-Schrodinger 
method will be proposed. We shall argue that an appropriate "optimal" normaliza- 
tion is equally well able to simplify many models containing N discontinuities in a way 
which enhances significantly the flexibility of the above-mentioned Morse-Feshbach 
single- well example. 

Our conjecture is based on several technical ingredients. Firstly, we imagine 
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that the contemporary computers shift the borderlines of feasibility of the initial 
zero order constructions. In this sense the standard solvable square well may easily 
be complemented not only by its textbook modification of finite depth [[| but by 
virtually any piece-wise constant potential V^°\r). This is discussed in Section 2. 
For definitness we pick up there the s— wave (£ = 0) rectangular or step-shaped 
example 

oo, r G (-oo,L )\J(L N+1 ,oo) 
Hj, re(L j ,L j+1 ), j = 0,l,...,N 

with impenetrable outer barriers H^+i = oo and H-\ = oo and with iV disconti- 
nuities and L > 0. These specifications are just convenient and could easily be 
altered. 

The study of perturbations of the more and more complicated solvable potentials 
([[]) reveals that the standard use of an unperturbed basis may become prohibitively 
cumbersome. The evaluation of all the necessary Rayleigh-Schrodinger overlap inte- 
grals does not seem rewarded by the resulting series 

ijj{ x ) = ^(o)( x ) + A + X 2 4j (2) {x) + ..., 

E = E<® + \E<» + \ 2 E& + . . . . 
A simpler recipe is asked for. In section 3 we propose, therefore, a new approach to 
corrections which circumvents the use of integrals. We shall see how it combines the 
continuity of our Schrodinger differential equation on certain finite intervals with an 
ease of their mutual matching. 

The formal appeal of our new technique lies in its unexpectedly coherent combi- 
nation of the matching of perturbed wave functions ip^ k '(r) with a "hidden" freedom 
of their normalization. The idea transcends its present application and makes the 
formalism quite universal. Our matching of perturbation corrections may be under- 
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stood as a useful alternative to the standard textbook recipe even in applications to 
smooth potentials. This is discussed in the last Section 4. 



2 Matching method in zero order 



We usually expect that a split of a given potential V(x) into a dominant part V^°'(x) 



and its perturbation A V^(x) simplifies our Schrodinger equation in its unperturbed 
limit A — > 0. Potentials V^°'(x) are predominantly chosen as harmonic oscillators. 
The more complicated shapes of V(x) can hardly be treated by perturbation expan- 
sions without resort to their discontinuous approximants. 

2.1 The piece-wise constant unperturbed potentials 

The first non-trivial s— wave step-like example (|l|) with N = 1, gauge H = 0, energy 
E = j3 2 , step Hi = j3 2 — 7 2 and abbreviations L = 0, L\ = P, L 2 = Q possesses the 
trivial wave functions 



Their matching at P fixes M and defines the energies as roots of the elementary 
trigonometric equation 



In the less common N = 2 example with L% = R let us abbreviate E = /3 2 , a 2 = 
Hi — (3 2 and 7 2 = (3 2 — H 2 and admit the complex a, (3 and 7 in 




x e (0, P) 




1 tg(3P = (3tg 1 {P-Q) . 
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In terms of B = (3P, C = j(Q — R), A = arctg(o;//3) and D = arctg(a/7) the 
spectrum follows from the similar matching condition 

e a(Q - p) cos(5 - A) cos(C + D) - e - a(Q - p) cos(5 + A) cos(C - D) = 0. (3) 

The graphical localization of its zeros is sampled in Figure 1 where we have chosen 
P — 1, Q = 2 and R = ir. With the double- well choice of H 2 = the Figure displays 
our double-well secular determinant at the four different heights of its central barrier. 
In principle, these curves range from zero up to the maximal k — (5 — y/Hi where 
they turn purely imaginary. Our picture shows just a small vicinity of the quasi- 
degenerate doublet of the two lowest energies. Due to the asymmetry of V^(x) 
their split only very weakly depends on the repulsive central core. 



2.2 Trigonometric symbolic manipulations 

After we move to the higher integers iV the assistance of a computer becomes wel- 
come. For example, the choice of iV = 4 may mimic a double tunneling. With 
H — H 2 — H 4 — and H 1 = H 3 = H > 0, at the energy E = k 2 = H cos a and 
with abbreviations L 4 = S, L 5 = T, 5 = y/H — E = \fH sin a > and 

F(a) = sin[(_R — Q) k(cx)] cos[(T — S) k(o) + a] cos[P k(o) — a], 

G(a) = sin[2o; — (R — Q) n(a)] cos[(T — S) k(cx) — a] cos[P k(o) — a] 
we derive the secular equation 



,(Q-P)S 



F{a)e( s - R ^ 5 + G(a)e 



-(S-R)S 



-(Q-P)s 



F(-a)e( s - R V + G(-a)e 



-(S-R)6 
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It is fairly transparent. In our last example with N = 6 the triplet of barriers in 

0, xe (0,P)\J(Q,R)\J(S,T)\J(U,W) 
V(QW) (x) = \ H, xe (P, Q) \J(R, S) U(T, U) 
oo, x e (-00, 0) u(W / , 00) 

requires similar strategy. With the four auxiliary functions 

F 1 =cos(P/t + a) sin[(P- Q) k + 2o\ sin[(T - S) k + 2 a] cos[(W - U) k + a], 
F 2 = cos{PK + a) sm[(R - Q) re + 2 a] sin[(T - S) re] cos[(W - U) re - a], 

F 3 = cos(Pre + a) sin[(P - Q) re] sin[(T - S) re] cos[(W - U) k + a], 
F 4 = cos(P k + a) sin[(P - Q) re] sin[(T -5)/t-2a] cos[( W - [/) re - a] 
of a and re = re(a) = re(— a) = VH cos a = y/E > and with the same 5 = \5\ -■ 
y/H — E = V^P sin a > as above, the exact secular equation is 

-Fi(a) e 2 ( p +*+ T ) 5 W + F 2 (a) e^+^+^W + F 3 (a) e 2i - p+s+T ^ 



-FA a) e ^+S+U)S(a) + F /_ a) e 2(Q+S+l/)*(«) _ Fj 



2 l-a e 



2(Q+S+T)8(a) 



-F 3 (-a) e 2 (Q+R+u)S(a) + F 4 (-oi) e 2(Q+H+r)«(a)l e [(--P-Q--R-s-T-i7)*(a)] =0 . 
The study of the other similar systems may be guided by this experience. The longer 
ansatzs remain tractable by computerized trigonometric manipulations. 



2.3 The general matching recipe 

Any Schrodinger bound state problem with a piece-wise constant potential is ex- 
actly solvable on each sub- interval (Lj, P,+i). Its two independent solutions are just 



the (trigonometric or hyperbolic) sines S • + ' ) (x) and cosines Cf } (x) determined in a 



»(+), 



unique way by the left initial conditions, 



C^\x) 
Si + \x) 



x=Li 



x=Lj 



= 1, 



0, 



d x c { ;\x) 



d x Sf\x) 



= 0, 



x=Lj 



(4) 



x=Li 
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They define the general solution simply as a superposition 

^°)(x) = cW(j)Cj +) (x) + rf( + )(j)5j +) (a;) , xe(L jt L j+1 ), j = 0, 1, . . . , N. (5) 

In the light of the obvious symmetry of our global problem on (L , L^ + i), one can 
equally well employ the alternative ansatz 

^W(x) = c (-)(j + l) cQ(x) + S-\ 3 + 1) S { r + \{ X ) , x G (L„ L i+1 ), (6) 
with the superscripted basis defined by the right initial conditions 
Cj+i i x ) 



= 1, d x C { r + \{x, 

x=Lj +1 



o. 



= 1. 

x=Lj +1 

We may immediately make the following two conclusions. 



0, d x S&(x 



x=Lj+i 



The standard matching of the logarithmic derivatives at every discontinuity 
x = Lj finds an elementary though rigorous guarantee in the obvious rule 

c (+) (j) = c (-) (j) = c(j)) d ( + ) {]) = d (-) (j) = d{j)i j = 0, 1, ... , N+l. 

This is the first simplification of our two alternative formulae (||) and ([]) for 
wave functions. 

In each interval (Lj,Lj + x), any two-point boundary conditions tp(°>(Lj) = c(j), 
ip^(Lj + i) = c(j + 1) make the solution unique. This statement is equivalent 
to the two simple rules 

c(j) = c(j + 1) Cj+\(Lj) + d(j + 1) Sj+\(Lj) , 
c(j + 1) = c(j) C\ + \L ]+1 ) + d{j) S\ + \L ]+l ) 
where, of course, c(0) = and c(N + 1) = and j = 0, 1, . . . , N. We may 
summarize that as a net result of our construction we have to solve just the 
2N + 2 homogeneous and linear algebraic equations for the 2N + 2 unknown 
parameters. 
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In a marginal comment let us also note that our choice of c(0) = and c(N + 1) = 
is dictated by the underlying physics which requires the most common Dirichlet 
form of the outer boundary conditions. In a more formal approach one could also 
keep in mind the possible alternative choices of the Neumann conditions d(0) = 
and d(N + 1) = or of the important periodic option for c(0) = c(N + 1) and 
d(0) = d(N + 1) etc. Unfortunately, their detailed study would already lead us far 
beyond our present task. 

2.4 Piece-wise polynomial unperturbed potentials 

We have seen that the practical use of our nonstandard matching recipe (§) is a well 
defined numerical problem. Its solution requires just the knowledge of the indepen- 
dent sine and cosine solutions and an evaluation of their values at all the points of 
the discontinuities x = Lj. It is important to note that our rigorous matching recipe 
works with the two independent bases. This enables us to avoid the more usual 
but, sometimes, less comfortable construction of the derivatives of the basis. In this 
sense, a slight generalization of this recipe may be easily applied to all the zero-order 
problems 

- -Sr ^ (0) (x) + V {0) (x) V (0) (x) = E<® ^ (0) far) (9) 
dx z 

with any piece-wise polynomial discontinuous potential 

p(±i) 

VV>{x±) ~ £ wf(j) (x ± - Lj) 1 , x + > Lj, x- < Lj. (10) 

1=0 

Locally (i.e., within the "double" intervals Jj = (Lj_ 1; L J+1 ) such that L_i = L 
and Ln+2 = L^+i) we may drop the redundant argument j and superscripts ^ and 
search for the exact wave functions in their respective left and right Taylor series 
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form 

M 

1>(a±) ~ ]T h±(j)(x+ ~ LjT, M - oo. 

ra=0 

In a purely numerical implementation Hodgson's tests [|l(| confirm the fast conver- 



gence of such a "local" recipe with N — > oo. In the present perturbation context 
the number of discontinuities x = Lj at j = (0), 1, 2, . . . , N, (N + 1) is fixed and, 
presumably, very small, N = 0(1). Still, in a compactified notation which parallels 
our previous p(±j) = construction we can drop the superscripts ± and re-write our 
local Taylor series as superpositions 

i/>(x) = c{j)C j (x)+d(j)S j (x) , xe(L j - 1 ,L j+1 ), j = 0,l,...,iV + l. (11) 

Their two components Cj(x) = C(x) and Sj(x) = S(x) are independent solutions of 
our ordinary differential Schrodinger equation @ again. They are uniquely deter- 
mined by their respective cosine-like and sine-like behaviour at x — Lj, 



C(x) 



Six) 



= 1, d x C(x) 

x=Lj 

= 0, d x S(x) 

x=Lj 



= o, 

x= Lj ^ 
= 1. 

x=Lj 



At the two outer boundaries L = L Q and i? = L^+i the current physical conditions 
ip{L) = ip(R) = acquire the most elementary form c(0) = and c(N + 1) = 
or, equivalents, c(l)&(L) + d{l) S^L) = and c{N)C N {R) +d(N)S N {R) = 0. 
This parallels again the above p(±j) = special case. The mutual matchings of the 
neighboring wave functions are all similar and we have the final compact set of the 
physical requirements 

c(j) Cj{L^ x ) + d(j) = c(j - 1), 

c(j) C 3 {L 3+l ) + d(j) S 3 (L J+1 ) = c(j + 1), (13) 
j = l,2,...,N. 
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We have to evaluate again the 4N input quantities Sj(Lj± 1 ) and Cj(Lj± 1 ) and solve 
the 2N— dimensional "secular" equation for the arbitrarily normalized coefficients in 
the local wave functions (|ll|) and for the global binding energy. All our piece-wise 



constant illustrative examples of subsections |2.1| and [E2] re-emerge after the choice 
of p(±j) = in potential (|T0| ) of course. 



3 Matching method for perturbations 

The separate 0(X k ) components of the perturbed Schrodinger equation have the well 
known non-homogeneous form 

- ^^ {k) (x) + V {Q \x) ^ {k) (x) -E {0) ^ k \x) = r {k - 1 \x) + E {k) ^ 0) (x) (14) 
ax 2 

with 

k-l 

T^ k - l \x) = -V {1 \x)^ k - 1 \x) + Y / E u) iJ (k ~ 1) (^), A; = 1,2,... . 

In principle, it defines the k— th corrections in terms of their predecessors ^^^'(x), 
ip( k ~ 2 \x), ...and E^ k ~ l \ E^ k ~ 2 \ "compressed" in the order-dependent right- 
hand side functions. 



3.1 Local solutions 



The matching method of subsection |2.3| does not use the (logarithmic) derivatives. 
This is rendered possible by a certain redundancy of our construction since domains 
Jj overlap. We shall now apply the same strategy to the implicit definition ([14]) of 
corrections E^ k > = e and ift( k '(x) = (p{e,£,x) split in four terms locally, 

^ k \ x ) = c «(j) cf\x) + S k \j) Sf\x) +eu(x) + ) ^\x), j = 1, 2, . . . , N. 

(15) 
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This is our key ansatz. Its two free parameters e and £ should facilitate the matching 
at the boundaries of J r Within each of these intervals and at e = and £ = the 
simplified order- dependent part of our non-homogeneous differential eq. (|14]) 

d 2 



dx 2 



+ V ( -°\x)-E^ 



(16) 



will define all the superpositions of the energy-independent functions C(x) and S(x) 
distinguished (and made unique) by the respective cosine- and sine-like initial condi- 
tions (Tl2|). The third auxiliary function u(x) = tOj{x) will be specified as a solution 
of a simpler, order-independent equation 

d 2 



dx' 



+ V(°\x)-E® 



uj(x) = ^ (0) (x) 



(17) 



with the different initial conditions 



UJj(x] 



x=Lj 



d x u)j(x] 



. 



The fourth component ip(°\x) is known. Its contribution is weighted by the last 
parameter £ chosen to shift the sum c^ k \j) + d^ k \j) by c^°\j) + d^°\j) at each j. 
Assuming that c^(j) + d^(j) ^ we may re-scale c^(j) + d^ k \j) = 1 at all k > 0. 

The k— and j— dependent variability of £ = £„■ does not violate the validity of 
our differential Schrodinger eq. ( |T3D locally. The global matching of its solutions flT5|) 
forms the last step towards an innovated perturbation method. 



3.2 Global solution in the k— th order 

Under the first nontrivial choice of N = 1 the matching of perturbation corrections 
degenerates to the left asymptotic-like boundary condition at x = L and its right 
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counterpart at x = R. This imposes the two linear algebraic constraints 

c<*> C{L) + (1 - c^) S(L) +eu{L) = 0, 
C W C{R) + (1 - c(*)) <S(i?) + e u(R) = 
upon the energy e = and coefficient c^(l) = 1 — d^^l) = X, 



S(L) 


-C(L), -u{L) 






S(L) 


S(R) 


-C(R), -uj(R) 


\ E ) 




S(R) 



One has to notice the possible absence of solutions of this system whenever its deter- 
minant vanishes. Such an apparent paradox just reflects an a priori open possiblity 
of degeneracy of the unperturbed spectrum. We only know a posteriori that the 
spectrum of the one-dimensional Sturm-Liouville problem cannot degenerate at all 

Ill- 
Let us abbreviate cf®(j) Cf ) {L i ) + S k \j) Sf\U) +eu}{L- l ) + & c^(i) = 

for N > 2 and, after the next choice of N = 2, contemplate the four independent 

boundary-and-matching conditions 

¥>i(Lo) = 0, <^(L 2 ) = <f 2 (L 2 ) [ = c( fc )(2)+e 2 c(°)(2)], 

^(Lx) = ^(Li) [ = c( fe )(l)+ac(°)(l)], <^ 2 (L 3 ) = 0. 
Both the re-normalization parameters enter these equations only in the form of their 
difference Z = Denoting c^(l) = 1-S k \l) = X and c^(2) = \-S- h \2) 

V we get the four linear relations 

XCf } (L ) + (1 -X)S?\Lo) +eu(L ) = 0, 
XCf ) (L 2 ) + (1 - X) 5? } (L 2 ) + su(L 2 ) + Zc^{2) = Y, 
YC { 2 k \L x ) + (1 - Y) <Sf(£i) H-ewCLO - Zc«>(l) = X, 
YC { 2 k \L 3 ) + (1 - F)5f(L 3 ) +^(L 3 ) = 
among the four unknowns e, X, Y and Z. This equation is easily solved by the 
four-by-four matrix inversion. 
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At an arbitrary N the general matching plus boundary formula comprises the 
2N equations 

X 3 C(L,_ 1 ) + (1 - X 3 ) S(Lj_i) + euj(Lj_i) = Xj-t + Zj c(°)(j - 1), 
Xj C(L j+1 ) + (1 - X 3 ) S(L j+1 ) + eu 3 (L 3+1 ) = X 3+1 - Z j+1 c^(j + 1), (19) 

j = l,2,...,N 

for 2N unknowns e, X 3 = c^ k \j) and Z 3+ i = ^ 3+ \ — £j. Our new perturbation 
prescription is complete. 

3.3 Illustration 

The detailed implementation of our matching formulae is straightforward. Its best 
illustration is provided by the solvable square well V^°\x) = V(sw)( x ) with the 
solvable constant perturbation V^'(x) = fi. In this extreme example the "survival 
of solvability" facilitates the understanding of formulae as well as a verification of 
their quantitative predictions without any use of a complicated algebra. For the sake 
of brevity we shall also pay attention to the N — 1 recipe in the first perturbation 
order only. 

3.3.1 Local wave functions 

In the first step it is easy to extract the particular solution = p(x) = 

(1/2) x cosx from the non-homogeneous differential eq. (0). Its order-dependent 
partner eq. ( JHty looks similar, 

d 2 



dx 2 
C(x 



+ V (0 \x)-E^ 



C(x) = t{x), 



= 1, d x C(x] 

x=X 



= 

x=X 
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S(x) 



S(x) = t(x 
d x S(x 



x=X 



x=X 



and possesses the similar particular solution —Qp(x). By means of the Ansatz 

oj(x) = P sinx + Q cosx + p(x), 
C(x) = U sinx + A cosx — Qp(x), 
S(x) = W sinx + B cosx — Qp(x) 

the initial conditions are easily satisfied by a suitable choice of the six optional 
constants P, Q, U, A, W and B. After an elementary trigonometry using the function 
q(x) = d x p(x) = (cosx — x sinx)/2 and an elementary orthogonal matrix 



R{x) 



COS X SIM 1 



— sin x cos x 



[R(-x)}- 1 



we get the result 





( 




= R(-x) 




V 



(p\ 

U j 

c + n P (x) 



R(-X) 



v -q{X) j 



\ 



( b\ 



( 



R{-X) 



\ W J 



Qp(X) 

\ d + n q (x) j 

Our local first-order solution is obtained by fully non-numerical means. 



3.3.2 Global matching and the energy 



Two-dimensional eq. ([18]) represents the physical boundary conditions at both ends 
of our interval of coordinates (L, R) . From its matrix elements 

u(L)=P, uj{R) = -(P + tt/2), C(L)=A, 
C(R) = VLit/2 - A, S(L) = B, S(R) = ttn/2-B 
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it is easy to deduce the answer = fi. This verifies the recipe and reproduces, 
incidentally, the exact result. 

It is instructive to notice that in the traditional Rayleigh-Schrodinger approach 
where the value of e is evaluated in advance our boundary conditions would be 
satisfied automatically In the present approach the variability of e is admitted 
breaking, in general, the boundary conditions 



. 



S(L)+euj(L) C(L)+eoo(L) I c« ^ 
S(R) +suj(R) C(R)+euj(R) \ d« y 

In an apparent paradox the latter equation seems quadratic (but is linear) in e. 
Immediate calculation reveals the above-mentioned energy correction uniquely, 



E (i) 



det 


S(L) C(L) 
S{R) C{R) 


fivr (B-A)/2 


det 


u(L) S(L) 
u(R) S(R) 


+ det 


C(L) oj{L) 
C(R) u(R) 


n(B-A)/2 



Wave functions may be reproduced in the similar manner. 



4 Summary 

We have seen in Section ||] that for many rectangular potentials V^°\x) the zero- 
order local wave functions are superpositions of elementary trigonometric functions 
sin f3x and cos fix with a real or purely imaginary argument. In conclusion we should 
now add that for an arbitrary piece-wise polynomial perturbation this reduces the 
construction of corrections in perturbation series (0) to an easy algebraic exercise. 
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4.1 Closed formulae for polynomial perturbations 



Firstly, let us notice that the inhomogeneous term r(x) in eq. ([H]) coincides with 
a certain superposition of products (x\k,l) = x k cos fix and (x\k, 2) = x k sin fix 
for polynomial perturbations. These functions may be denoted and treated as a 
partitioned basis { \k,j) }y=i,2 with fe = 0, 1, . . .. In this basis the action of the 
unperturbed differential operator 

d 2 



n 



+ V^\x)-E^ 



has a closed explicit form 



(2/3) 1 H x cos f3x = sin fix, 
(4/3 2 ) _1 7i (/3x 2 cos /3a; — x sin /fa) = x sin /fa, 
(12/3 3 )- 1 ?? (2fi 2 x 3 cos /fa - 3/5 x 2 sin /fa - 3x cos fix) = x 2 sin fix, 
(8/3 4 ) _1 7i (/3fa 4 cos /fa — 2fi 2 x 3 sin /fa — 3fi x 2 cos /fa + 3 x sin /fa) = x 3 sin /fa, 



and, mutatis mutandis, for cosines. Thus, any practical computation will immedi- 
ately generalize our previous trivial constant-perturbation example. 

In the constructive proof of the latter relations we firstly represent the action of 
our operator 7i on each element of the basis as a superposition of the other basis 
states. It is easily shown by explicit differentiation that the coefficients of these 
superpositions form an infinite matrix Q with the mere three nonzero diagonals. In 
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the second step, we introduce a two-by-two partitioning of the matrix Q and denote 



Q 



\ 



o ... 

b x ... 

c 2 b 2 ... 

c 3 6 3 ... 

c 4 W ... 

V J 

The submatrices bk = 2k/3a and = —k(k—l)I are elementary and two-dimensional, 





r 


1^ 




fi 




(= v^T) = 














v- 1 








1 J 



In the third step we verify that the /e/£ inverse Q L of our singular matrix still exists 
and has the elementary form 

/ 



Q 1 



b^ 1 

-b 2 l c 2 b^ 







&2 1 








b^c z b 2 l c 2 bl l ~b^c 3 b 2 l 63 1 



V : : : : '• 7 

The explicit form of its two-by-two submatrices remains fairly compact, 



( st )„ +2 „ +1 = < 26 >~ 3 <" + 2 > ff - ( sL )„ +3 „ +1 = < 26 »" 4 (" + 2 ><" + 3 > 7 

(fi 1 ) ^ x, = -(2*)" 5 (n + 2)(n + 3)(n + 4)<r, . . . , n = 0, 1, . . . . 

\ / n+in+1 

This completes the proof. The separate rows of the non-partitioned matrix Q L 
determine the particular solutions of our fundamental differential equations (fL6"D and 

(ED- 
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4.2 Generalizations 

Our last observation was extremely pleasant and encouraging. It immediately im- 
plies that for the piece-wise constant unperturbed potentials our new perturbation 
construction remains non-numerical for each polynomial perturbation. In this sense 
the usual start from a harmonic oscillator may find here an unexpectedly feasible me- 
thodical alternative even in analyses of continuous models. We have seen that up to 
the discontinuities at the lattice points x = Lj, j = 1, 2, . . . , N the one-dimensional 
or s— wave functions remained basically non-numerical. 

Technical complications may emerge beyond s— waves, for polynomial Vj (x) 
and for the nonpolynomial perturbations. All these problems may appear quite 
naturally in many applications. In such a case both the unperturbed problem and 
the evaluation of corrections become much more numerical. Still, an implementation 
of our perturbation recipe remains virtually unchanged, consisting of the following 
six steps. 

• SI. We solve the unperturbed differential Schrodinger equation with the ap- 
propriate initial conditions QT2] ) in all the domains Jj. Their number N is a 
fixed and, presumably, very small integer parameter. 

• S 2. We solve the linear algebraic system of the 2 N homogeneous equations 
(|i"5D. This determines the unperturbed energy as well as the unperturbed 
norms cf^ and dj of the matched zero order wave functions. 

• S3. We solve all the auxiliary initial- value problems and generate the 
N functions u(x) = Uj(x). In particular, their values tjj(Lj±i) have to be 
computed at the boundaries of all domains. 
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• S 4. In the given order k = 1,2,... and in every domain Jj we specify the 
"input" finite sum r^^ix) and solve the doublet of the initial value problems 
fll6|). This determines the functions Cf\x) andsf (x) as well as their special 
values Cf\L j±l ), Sf\L j±l ). 

• S 5. We solve, finally, our finite set of the 2N "effective" or "model-space" 
linear algebraic eqs. (p~Q|) . This defines the k— th energy correction e = 

the N matched norms c^ k \j) and the N — 1 local re-normalization parameters 
Zj = ^ — Q^i in the wave functions. 

• S 6. If needed, we move to the next order k and return to step S 4. 

We may conclude that in the future applications of our new perturbation prescrip- 
tion its present algebraic (trigonometric) exemplification may be complemented by 
some local versions of the current semi-analytic Taylor series constructions or by the 
various discrete (e.g., Runge Kutta) purely numerical implementations etc. 

Acknowledgements 

Partially supported by the grant Nr. A 1048004 of the Grant Agency of the Academy 
of Sciences of the Czech Republic. 

Figure captions 

Figure 1. Momentum dependence of secular determinant @) for four different barriers 
H x = 10, 15, 20 and 25. 
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